Transport through a double-quantum-dot system with noncoUinearly polarized leads 
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We investigate linear and nonlinear transport in a double quantum dot system weakly coupled 
to spin-polarized leads. In the linear regime, the conductance as well as the nonequilibrium spin 
accumulation are evaluated in analytic form. The conductance as a function of the gate voltage 
exhibits four peaks of different height, with mirror symmetry with respect to the charge neutrality 
point. As the polarization angle is varied, due to exchange effects, the position and shape of the 
peaks change in a characteristic way which preserves the electron-hole symmetry of the problem. 
In the nonlinear regime various spin-blockade effects are observed. Moreover, negative differential 
conductance features occur for noncoUinear magnetizations of the leads. In the considered sequential 
tunneling limit, the tunneling magnetoresistance (TMR) is always positive with a characteristic gate 
voltage dependence for noncoUinear magnetization. If a magnetic field is added to the system, the 
TMR can become negative. 

PACS numbers: 72.25.-b, 73.23.Hk, 85.75.-d 
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I. INTRODUCTION 



Spin-polarized transport through nanostructures is at- 
tracting increasing interest due to its potential applica- 
tion in spintronicsii^ as well as in quantum computing^. 
Downscaling magnetoelectronics devices to the nanoscale 
implies that Coulomb interaction effects become increas- 
ingly importanli^. In particular, the interplay between 
spin-polarization and Coulomb blockade can give rise to 
a complex transport behavior in which both the spin 
and the charge of the "information carrying" electron 
play a role. This has been widely demonstrated by 
many experimental studies on single-electron transis- 
tors (SETs) with ferromagnetic leads, with central ele- 
ment being either a ferromagnetic particle^'^'*, normal 
metal particle o^'^" , a two-level artificial moleculei"^, a 
Ceo moleculei^, or a carbon nanotubei^, showing the 
increasing complexity and variety of the investigated 
systems. Initially, the theoretical work was mainly fo- 
cused on the difference in the transport properties for 
parallel or antiparallel magnetizations in generic spin- 
valve SETsiiOSiieaiiisaSiao^^^^, More recently, 

the interplay between spin and interaction effects for 
noncoUinear magnetization configurations has attracted 
quite some interest both in systems with a continuous 
energy spectru m^^i^^i^^'^^ , as well as in single-level quan- 
tum dots2£i^i2^i^i^i^i2&i^, many-level nanomagnets^ 
and in carbon nanotubes quantum dots^S. In the non- 
coUinear case, a much richer physics is expected than 
in the coUinear one. For example, two separate ex- 
change effects have to be taken into account. On the 
one hand, there is the non-local interface exchange, in 
scattering theory for non-interacting systems described 
by the imaginary part of the spin-mixing conductance-^, 
and which in the context of current-induced magnetiza- 
tion dynamics acts as an effective field^"'^. Such an ef- 
fective field has been found experimentally to strongly 
affect the transport dynamics in spin valves with MgO 
tunnel junctionSi^. This effect has recently also been in- 




FIG. 1: (Color online) Schematic picture of the model: A 
double quantum dot system attached to polarized leads. The 
significance of the on-site and inter-site interactions U and 
V, respectively, is depicted. The source and drain contacts 
are polarized and the direction of the magnetizations is 
indicated by the arrows. 



volved to explain negative tunneling magnetoresistance 
effects in carbon nanotube spin valves^^ and called spin- 
dependent interface phase shifts^^'^'^. The second ex- 
change term is an interaction-dependent exchange effect 
due to virtual tunneling processes that is absent in non- 
interacting systems^^'^^''*-:-- . This latter exchange effect 
is potentially attractive for quantum information process- 
ing, since it allows to switch on and off magnetic fields 
in arbitrary directions just by a gate electric potential. 

Recently, there has been increasing interest in dou- 
ble quantum dot systems realized e.g. in semiconduct- 
ing structures*^ or carbon nanotubes^, as tunable sys- 
tems attractive for studying fundamental spin correla- 
tions. In fact the exchange Coulomb interaction in- 
duces a singlet-triplet splitting which can be used to 
perform logic gates^. Moreover, Coulomb interaction 
together with the Pauli principle can be used to in- 
duce spin-blockade when the two electrons have triplet 
correlations^2i^i^i^. The Pauli spin-blockade effect can 
be used to obtain a spin-polarized current even in the ab- 
sence of spin-polarized leads; it requires a strong asym- 
metry between the two on-site energies of the left and 
right dot. 

So far transport through a DD system with 
spin-polarized leads has been addressed in few 
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theoretica P^i^^'^^ and experimental^^ works, for the case 
of collinearly polarized leads only. While Rei2^ addresses 
additional Pauli spin-blockade regimes when one lead is 
half-metallic and one is non-magnetic, Ref;^'' focusses on 
the effects of higher order processes in symmetric DD 
systems, which can e.g. yield a zero bias anomaly or 
a negative tunneling magnetoresistance. Inii Coulomb 
blockade spectroscopy is used to measure the energy dif- 
ference between symmetric and antisymmetric molecular 
states, and to determine the spin of the transferred elec- 
tron. 

In this work we investigate spin-dependent transport 
in the so-far unexplored case of a double-dot (DD) sys- 
tem connected to leads with arbitrary polarization di- 
rection. Specifically, we focus on the low transparency 
regime, where a weak coupling between the DD and the 
leads is assumed. Our model takes into account interface 
reflections as well as exchange effects due to the interac- 
tions and relevant for noncoUinear polarization. We focus 
on the case of a symmetric DD, so that rectification ef- 
fects induced by Pauli spin-blockade are excluded. In the 
linear transport regime the conductance is calculated in 
closed analytic form. This yields four distinct resonant 
tunneling regimes, but due to the electron-hole symmetry 
of the DD Hamiltonian, each posses a symmetric mirror 
with respect to the charge neutrality point. However, 
by applying an external magnetic field, this symmetry is 
broken, which can lead to negative tunneling magnetore- 
sistance features. Finally, in the nonlinear regime some 
excitation lines can be suppressed for specific polariza- 
tion angles, and negative differential features also occur. 

The method developed in this work to investigate 
charge and spin transport is based on the Liouville equa- 
tion for the reduced density matrix (RDM) in lowest or- 
der in the reflection and tunneling Hamiltonians. The 
obtained equations of motion are fully equivalent to those 
that could be obtained by using the Green's function 
method'^^:^^ in the same weak-tunneling limit. The ad- 
vantage of our approach is that it is, in our opinion, eas- 
ier to understand and to apply for newcomers, as it is 
based on standard perturbation theory and does not re- 
quire knowledge of the nonequilibrium Green's function 
formalism. 

The paper is organized as follows. In Sec. [IT] we intro- 
duce the model system for the ferromagnetic DD single- 
electron transistor. In Sec. IIIII the coupled equations of 
motion for the elements of the DD reduced density ma- 
trix are derived. Readers not interested in the derivation 
of the dynamical equations can directly go to Sees. IIVI 
andlVl where results for charge and spin transfer in the 
linear and nonlinear regime, respectively, are discussed. 
Finally, we present results for the transport characteris- 
tics in the presence of an external magnetic field in Sec. 
IVIl Conclusions are drawn in Sec. IVIII 



II. THE MODEL 

We consider a two- level double-dot (DD), or a single 
molecule with two localized atomic orbitals, attached to 
ferromagnetic source and drain contacts and with a ca- 
pacitive coupling to a lateral gate electrode. The system 
is described by the total Hamiltonian 

H = Hq+Hs+HcI + Ht + Hr., (1) 

accounting for the DD Hamiltonian, the source (s) and 
drain (d) leads, and the tunneling and reflection Hamil- 
tonians, respectively. The two contacts are considered 
to be magnetized along an arbitrary, but fixed direction 
determined by the magnetization vectors rfia. The two 
magnetization axes enclose an angle Q G [0°,180°] (see 
figure [T]). The spin quantization axis Za in lead a is par- 
allel to the magnetization rfia of the lead. The majority 
of electrons in each contact will then be in the spin-up 
state. The Hamiltonians Hs,Hd that model the source 
(s) and drain (d) contacts read {a = s, d) 

Ha = J2^£ka^ - Ma)4fc^^C„;,^^, (2) 

where cj^^,^ and c^^^ are electronic lead operators. 
They create, respectively annihilate, electrons with mo- 
mentum k and spin CTq in lead a. The electrochemical 
potentials fia — fJ-Oa + ^Va contain the bias voltages Vs 
and Vd at the left and right lead with Vg — Vd = Vhias- 
There is no voltage drop within the DD. We denote in 

the following £fecr„ - Ha ■= £aka^- 

Tunneling processes into and out of the DD are de- 
scribed by Ht- We denote with rf^cr^j '^ao-c ^he creation 
and destruction operators in the DD. We assume that 
tunneling only can happen between a contact and the 
closest dot, so that we can use the convention that to the 
leads indices a = s,d correspond a = 1,2 for the DD. 
With ta the tunneling amplitude we find 

Ht=Y1 (^"4<T„c„fc<,^ -f t^ci^^jaaj- (3) 

The so-called reflection-Hamiltonian includes re- 
fiection events at the lead- molecule- interfac o^^'^^ . For 
strongly shielded leads the overall effect is the occurrence 
of a small energy shift A/j, induced by the magnetic field 
in the contacts and built up during several cycles of re- 
flections at the boundaries. It reads 

Hn = -An (dlu^U " <lJaJ- (4) 

a—s,d 

Finally, the DD Hamiltonian needs to be specified. As 
spin quantization axis of the DD, zq, we choose the direc- 
tion perpendicular to the plane spanned by Zs and Zdr^ 
(see Fig. [5]). The two remaining basis vectors ^0 and ijQ 
are along Zg + Zd, respectively z^ — Zd- The matrices which 
mathematically describe the above transformations read 
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FIG. 2: (Color online) The spin quantization axis of the 
double-dot, zq, is chosen to be perpendicular to the plane 
spanned by the magnetization directions rhs, fhd in the leads. 
The latter enclose an angle O. 



normalized to one, the condition ctg + /?q = 1 holds. The 
energy difference between the S = groundstate and the 
triplet is given by the exchange energy 

J =]^{/^-U + V) = 2|6|(i?+ Vl + i?') , 

where A = A\h\VlTW and R = (U - V")/(4|6|). 
Besides the triplet, one observes the presence of higher 
two-particles excited states with total spin 5 = 0. 



Abbr. State Eigenvalue Spin 



1 Z+e+^e/* +e-*Q/4 
71 1^ _e+»e/4 +e-»e/4 



M. 



(5) 



We express the DD Hamiltonian in the localized ba- 
sis, such that e.g. |-|-, — ) describes a state with a spin-up 
electron on site 1 and a spin down electron on site 2 (with 
spin directions expressed in the spin-coordinate system of 
the DD). Such state can be obtained by applying creation 
operators on the vacuum state, i.e., ) = (i^(i2j,|0). 
In general, the ordering of the creation operators is de- 
fined as (i^(i2|(ijj^(i2j0). The DD Hamiltonian then 
reads 



(6) 



^0 = ^'^cd'l^^d^^^^+b^idl^^d^^^+dl^^di^ 

Q(T0 CT0 

a— 1 (Tq 

2 



where the spin-index indicates that the operators are 
expressed in the spin-coordinate system of the DD. The 
tunneling coupling between the two sites is 5, while U 
and V are on-site and inter-site Coulomb interactions. In 
the remaining we consider a symmetric DD with equal 
on site energies Si = £2- Thus we can incorporate the 
on-site energies in the parameter ^ proportional to the 
applied gate voltage Vgate- 

To understand transport properties of the two-site sys- 
tem in the weak tunneling regime, we have to analyze 
the eigenstates of the isolated interacting system. These 
states, expressed in terms of the localized states, and the 
corresponding eigenvalues are hsted in table The 
table also indicates the eigenvalues of the total spin op- 
erator. The groundstates of the DD with odd particle 
number are spin degenerate. In contrast, the ground- 
states with even particle number have total spin S — 
and are not degenerate. In the case of the two particles 
groundstate the parameters ag determine whether 

the electrons prefer two pair in the same dot or are delo- 
calized over the DD structure. Since the eigenstates are 
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A = 4\b\VT+W, 
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TABLE I: Eigenstates of the double-dot system and corre- 
sponding eigenvalues and parity. In the limit \b\ —> 00, where 
the inter-dot hopping is unhindered, 7? — > and ao —> Po- 
For l&l — > 0, i.e. no inter-dot hopping takes place, we find, if 
U > V, that R +00 and qq ^ 1, /3o ^ 0; the state |2) 
then becomes degenerate to j2'(0)), forming a Heitler-London 
state. In turn, if U < V then R —00 and ao —> 0, l3o ^ 1. 

Finally, we remark that and Hj^ contain operators 
of the DD, and d^^^, with spin quantization axis 



4 



of the leads, while Hq is already expressed in terms of 
DD operators d]^^^ and d^^^ with spin expressed in the 
coordinate system of the DD. 



III. DYNAMICAL EQUATIONS FOR THE 
REDUCED DENSITY MATRIX 

In this section we shortly outline how to derive the 
equation of motion for the reduced density matrix (RDM) 
to lowest non vanishing order in the tunneling and re- 
flection Hamiltonians. The method is based on the well 
known Liouville equation for the total density matrix in 
lowest order in the tunneling and reflection Hamiltonian. 
Equations of motion for the reduced density matrix are 
obtained upon performing the trace over the leads de- 
grees of freedom^'^, yielding, after standard approxima- 
tions, Eqs. (dni) and below. In the case of spin- 
polarized leads, however, it is convenient to express the 
equations of motion for the RDM in the basis which di- 
agonalizes the isolated system's Hamiltonian and in the 
system's spin quantization axis. After rotation from the 
leads quantization axis to the DD one, Eq. I|2ip. which 
forms the basis of all the subsequent analysis, is obtained. 

Let us start from the Liouville-equation for the total 
density matrix {t) in the interaction-picture 

,n'^^[H^^{t) + Hi,{t),p\t)l (7) 

with Ht and transformed into 

the interaction picture by H;!p^j^{t) = 

where Iq indicates the time at which the perturbation 
is switched on. Integrating ([7]) over time and inserting 
the obtained expression in the r.h.s. of ([7]) one finds 
equivalently 

P'it) = ~^[Hkit),p'ito)] - ^mit),p'ito)] (8) 

The time evolution of the reduced density matrix (RDM) 

pl{t):=Tri,adsip\t)) (9) 

is now formally obtained from ([5]) by tracing out the lead 
degrees of freedom. To proceed, we make the following 
standard approximations: 

i) The leads are considered as reservoirs of non- 
interacting electrons which stay in thermal equilibrium 
at all times. In fact, we only consider weak tunneling 



and therefore the influence of the DD on the leads is 
marginal. Hence we can factorize the density matrix of 
the total system approximatively as 

p'it) « p^(t)p,pd, ^ (10) 

where ps and pd are time independent and given by the 
usual thermal equilibrium expression for the contacts 

Ps/d = z d — ' ^'^^^ P being the inverse temperature 
and ^s/d the partition sums over all states of lead s j d. 

ii) We consider the lowest non vanishing order in Ht/b - 

iii) We apply the Markov approximation, i.e., in the 
integral in Eq. ([5]) we replace PqII') with PQ{t). In other 
words, it is assumed that the system looses all memory 
of its past due to the interaction with the leads electrons. 

Furthermore, being interested in the long term behav- 
ior of the system only, we send to — > — oo. We finally 
obtain the generalized master equation (GME) for the 
reduced density matrix 

PoW = ~^Tri,ads[Hj,it),pl,{t)p,p4 (11) 

- ^J^ dt''Tneads{[HUt),[Hkt-n,PUt)PsPd]])- 

A. Contribution from the tunneling Hamiltonian 

In the following, we derive the explicit expression for 
the GME in the basis of the isolated DD. For simplicity 
we omit the contribution of the reflection Hamiltonian 
in a first instance. When we shall have obtained the 
final form of the GME due to the tunneling term, we 
will see that it is easy to insert the contribution from 
the reflection Hamiltonian. Let us then start from the 
tunneling Hamiltonian in the interaction picture 

Hkt) = E E (12) 

akaa i.j 

|i)(j|exp [i{ )t/n] + h.c. 

where {d^aj,^ = (Adaajj) and (^LJ^^. = (iMLji) 
are the electron annihilation- and creation-operators in 
the spin-quantization axis of lead a expressed in the ba- 
sis of the energy eigenstates of the quantum dot sys- 
tem. To simplify (|lip standards approximations are in- 
voked, i) The first one is the secular approximation: 
Fast oscillations in time average out in the stationary 
limit we are interested in, and thus can be neglected. 
Together with the relation TrieadsiPsPdclka^'^a'k'a'J = 
Skk'Saa'Saa'faieaka), whcrc faisakcr) IS the Fcrmi func- 
tion, and the cyclic properties of the trace we get 
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^oW = -Jp / dfY, \t^\'{ (13) 

+ ^ /a(eQfeCT„)(dQa„)f ('^la„);mN)("^|P0(i) SXp [i(em - £; + eak<y„)t" /K\ 
ilm 

+ ^ Ia{^aka^)){d'^aaJtMaaJlrn\i){'m\PQ)it) ^Xp [-«(£; " Em + £akcrjt"/h] 

ilm 

+ ^fa{£akcr^)PQ{t){d^^Ju{dl,^Jlm\i){m\ex.p [-iiSi ~ £l + £aka^)t"/h] 
ilm 

+ ^ f<^i'^'^kaJ)pQ{t){d'l„Jii{d^^Jijn\i){m\exp [+i{ei - Si + £aka^)t" /h] 

ilm 

- ^(1 - fa{£aka^)){daaJi]PQ{t)]l{d'^aaJl7n\i){'rri\eyiY> [+i(e„i - £l + £aka^)t" / h] 
iljrn 

- X! •/'"('^"*;<T„)((i]j^^)y/5Q(t)j/(d„^^);„i|i)(m|exp - £,„ + £afc<T„ )i'7^] 

- ^ /a(eafe^„))Ko-JyP0Wji(4o-JimN>(»T^|exp [-i(£j - £j + Eckajt" 
iljra 

- XI /"(Sofe<T„)(dL„)»jPoWi/Ka„)imN>(»7i|exp [+i(£j - £» + £afe<T„ )i'7^] }• 

I 



ii) For the second approximation we notice that we wish 
to evaluate single components {n\pQ\m) of the RDM in 
the system's energy eigenbasis. Therefore, we assume 
that the DD is in a pure charge state with a certain num- 
ber of electrons TV and energy En ■ In fact, in the weak 
tunneling limit the time between two tunneling events is 
longer than the time where relaxation processes happen. 

I 

(a) 
(b) 
(c) 

(d) 

(e) 
(f) 

In dH]) we used the notation (n|pQ^"^|m). 



That is, we can neglect matrix elements between states 
with different number of electrons, and only regard el- 
ements of Pq which connect states with same electron 
number N and same energy En- So we can divide /3q 
into sub-matrices labelled with N and En and find 



I 

By convention, {X); z'' Sj' Sft, /t'} means that in each 
line (a)-(f) we sum over the indices occurring in this 



■EmN 

rnm 



E 



l,l'i£\N-l) je\ENN) h,h'e\N+l) 



TT J £k-£h+£j 



fai£k)Daa^{ek) 



(dca^ )nh (dL„ )h3Pf^^{t) 



(1 - fa{£] - £/)) Daa^ {£j " £/) 



fa{£h - £j)Daa^ {£h " £j) ~ ~ j 



(1 - fa{£] - £/)) Daa^ {£ j - £/) 



dei 



(1 - fa{£k)) DaaA£k) 



idLjnlid^.Jljpf^''it) 



fa{£k)Daa^ (£/c) 



£k~ £h + £j 

(1 - fa{£k)) Dacr^{£k) 



nj' ii)idaa-^)jhidli^^^)hm 



d£k 



£k - £j + £l 



Pnj 



-2 (1 - fa{£h - £j)) Daa^ {£h - £j)ido,aJnh'{dl^JhmPh'\^^\t) 

-2Ue, - £i)D^aA^, - £i){di„Jnv{d^.Ji^pf,f-\t) } . 



(14) 
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line only. Notice that the sum over j is restricted to 
states of energy Ej = — En = E^- For the 

states with iV ± 1 electrons, we have to sum over all 
energies, therefore we indexed the density matrix with 
Eh = Eh' respectively Ei = Ei' in lines (e) and (f). 
Further, we replaced the sum over k by an integral: 

J2k * / dSaka„Daa^{eakcr^), wherC Daa^{£aka„) de- 
notes the density of states in lead a for the spin direction 
CTq,, and applied the useful formula 



For later convenience we also define 



/' G 
dSakcr 7 



-E)t" 



G{£akij^) 



(15) 



where the prime at the integral denotes Cauchy's 
principal part integration. In our case G{eakaa) ~ 

Daa^i£akcr„)fai£akaJ with /+ = fa and / " = 1 - . 

In order to simplify the notations we replaced Eakaa, by 
Sk in (m. 



B. Transformation into the spin coordinate system 
of the double-dot 



In the previous section we introduced the transforma- 
tion rules for changing from the lead spin coordinates CTq, 
into the DD spin coordinates CTq. These rules give 




1 



V2\ 



+e 



"To 



(16) 

with :=-f, ed-=+^- 

Thus, Eq. p4)) can be easily expressed in the DD spin 
quantization axis. For example it holds 



(17) 



2 I Da+^ /± {E) - Da-Jt [E] aQ + 4 



0' 



and its related principal part integral 



Pt^o'^^E) := fdeF^^^^^,Je){e~E)-\ 
Contribution from the reflection Hamiltonian 



In order to give the full expression for the GME in the 
system's eigenbasis, we need to compute the contribution 
from the reflection Hamiltonian in Eq. (jlip . In analogy 
to what we did to evaluate the contribution from the tun- 
neling Hamiltonian, we must first transform Hfi into the 
interaction picture and then perform the secular approx- 
imation to get rid of the time-dependence. To start, we 
express Hf( in the DD spin quantization basis. 



Hi 



ie\N-i) 



(18) 



The commutator is easily evaluated to be 



--Trieads[HR, PQ{t)psPd\ = 

'^E^« E E E *-o-o 



(19) 



[dl 



d.',Jj){j\pUt)-pUt)di,^,d^',,;j){j\ 



In order to include this commutator in the master equa- 
tion let us introduce the abbreviation 



where we introduced 



Now we can add RaaQu'^ in (fT4| in the lines (6) and {d) 
to find the final form of the complete master equation in 
the DD spin-coordinate system. It reads 



= E E E ' E ' E { (^d 

a=sA (T0,crO' [(.('elW-l) je\ENN) hJi'e\N+l) J 
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(b) 



(e) 
(f) 



F , (si - ei) P , (si - ei) + Ra 



QCTq 



' hm 



5^0 '^O 



70 ^j-;' 



-2$ 
-2$ 



hm 



, F+ 



D. The current formula 



We observe now that ([2T|l can be recast in the Bloch- 
Redfield form 



We solve Eq. (|22p numericaUy and use the result to eval- 
uate the current flowing through the DD, as shown in 
the forthcoming sections. At low bias voltages, however, 
we can make some further approximations to arrive at an 
analytical formula for the static DC current. 



hh' 



(22) 



where the sums in (|22p run over states with fixed par- 
ticle number: £ {\EnN)}, h,h' e {|A^+1)}, 
1,1' £ {|iV— 1)}. The Redfield tensors are given by 



j^NN - ST \^ \x {t.{+)NN+1 U+)NN-1 
^nrnjj' — ["mj' yi^ a,nhhj ~^ a,nUj 



a l(or h) 

"'J V a,j' nnrn a.j' lim 



T^NN±l ^ \p /-^{+)NNt1 , p(-)AfJVTl 
nmkk' / j \ a,k'7nnk a.k' 



ik y ' 



(23) 



(24) 



where the quantities ^l^njjk^^ easily read out 

from (HH). They are 



-.(±)AfA'+l 

a,nhh'k / , i- ^ - ""©"q i i L-a(T0Crg' 

i 



'■©<7© 



With the stationary density matrix p^gi being known, 
the current (through lead a = s/d — ±) follows from 

A'" n,n',j 

(25) 



IV. THE LOW-BIAS REGIME 

A. General considerations 

A low bias voltage ensures that merely one channel is 
involved with respect to transport properties. Here we 
focus on gate voltages which align charge states N and 
iV-|-l. Moreover, we can focus on density matrix elements 
which involve the energy ground states E^^^ and 
only. In the following we shall use the compact notations 



'N 

N " 



4N) 



rnjn 



(26) 



Evaluation of the current requires the knowledge of 
/5^^ and p'"q'^^\ i.e. a solution of the set of coupled 
equations obtained from (pij) . or, equivalently, from (|22p . 
In the low bias regime this task is simplified since i) terms 
which try to couple states with particle numbers unlike 
N and N + 1 can be neglected; ii) we can reduce the 
sums over /i, h' in the equation for /5|^\ and over 1,1' in 
the equation for P^q'^^'' to energy-groundstates i?^^ and 

-Sj^^j^, because all the other transitions are suppressed 
exponentially by the Fermi-function. Notice, however, 
that these two approximations are not appropriate for the 
principal-part-terms, since they are not energy conserv- 
ing. The resulting equations for p[^^ and p[^~^^\ Eqs. 



(jBl[) and (jB2p respectively, can be found in appendix IB] 
In the following we shall apply those equations to derive 
an analytical expression for the conductance in the four 
different resonant charge state regimes possible in a DD 
system, i.e.. 



N = 0^ N = 1, 
N = 2 ^ TV = 3, 



N = 1^ N ^2, 

iV = 3 ^ TV = 4. (27) 



In all of the four cases we get a system of five coupled 
equations involving diagonal and off-diagonal elements 
of the RDM. The matrix elements of the dot operators 
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between the involved states entering these equations are 
given in appendix [Xj Before going into the details of 
these equations, it is instructive to analyze the structure 
and the physical significance of the involved RDM ele- 
ments. 



B. The elements of the reduced density matrix 

N = 0. In the case of an empty system we have only 
one density matrix element in the corresponding block 
with fixed particle number A'^ = 0, i.e., 



(28) 



describing the probability to find an empty double-dot 
system. 

N = 1. In this case we have four eigenstates for the 
system, where the two even ones build the degenerate 
groundstate and the two odd ones are excited states (see 
tabled}. In the low-bias-regime we only need to take into 
account transitions between groundstates. Therefore we 
have to deal with the 2 by 2 matrix 



(29) 



\^ wie-'"'^ Wii J 

The total occupation probability for one electron is 

Wi ■.= Wi-^ +Wii. (30) 

The meaning of the off-diagonal elements, the so called 
coherences, becomes clear if we regard the average spin 
in the system 

5f) = iTr(ar"'VL'^(t)) , (31) 

where i = x,y,z and cr^^""'* are the Pauli spin-matrices. 
This yields 



Wicosofi, S'^"'^-' = — Wi sinai, (32) 



(33) 



N = 2. For the case N = 2we actually have six different 
eigenstates, but only one of them, |2), is a groundstate 
(with spin 5 = 0), see table H] Only this groundstate 
must be considered in the low-bias-regime, yielding 



(34) 



This element describes the probability to find a dot with 
two electrons. 

N = 3. In this case we have again four eigenstates for 
the system, whereas the two odd ones build the degen- 
erate groundstate and the two even ones are excited. In 



the low-bias-regime we only need to deal with the 2 by 2 
matrix involving the three-particle groundstates 

The total occupation probability for three electrons is 

W3 ■.= W3^ + W3i. 

As for the case N — 1, the off-diagonal elements yield in- 
formation on the average spin S^^^ — \Tr ^crf ""'V©^ (O) 
in the system through the relations 

S''^ = =W3COsa3, S^^^ — —wssinas, (36) 

(37) 



N = 4. Finally, if the double quantum dot is com- 
pletely filled with four electrons we only have one non- 
degenerate state. Correspondingly, there is only one rel- 
evant RDM matrix element. 



(38) 



describing the probability to find four electrons in the 
system. The total spin is S = 0. 

Hence, we see that in all of the four cases ([77)) we get a 
system of five equations with the five independent phys- 

(i) (i) (i) 

ical quantities Wn,Wn+i and Sx \ Sy \ Sz with i = 1 
or 3. 



C. The conductance formula 

We shall exemplarily present results for the resonant 
transition = 1 = 2. For the other transitions 

similar considerations apply. The quantity of interest are 
Wi,W2,si^\sP,si'^\ related through Eqs. 1^, JMl) 
to the density matrix elements of /5q\ From Eqs. (|Bip . 
(|B2[) . and the table Hill of the appendix we finally obtain 
with Wi = l- W2, 



Wi 



a—s,d 

2F+,,{p2)W,-4F-,,{p2)W2 
4F+^{p2)S^'^-m^) , 



(39) 



|2;„2 



2F+^{P2)S^^'> 



[f+ ^{p2)W, - 2F-^{p2)W2)m^ 



Trfc 



(40) 



We have introduced the notation 4fc^ :— (ao ± /3o)^- AH 
the non-vanishing principle-value factors Pa-\[ and the 
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reflection-parameter Ra]i have been merged to the com- 
pact form 



Vo^ifii, {E2} - e[°^) := --[P- ^(^1) + R^^^] 



(41) 



-Pq^i (^2" — £le) — fc^-Pii (£2"' — £le) , 



where we introduced the chemical potential fiN+i 



E' 



(0) 
JV+l 



i?^^ and {E2} denotes the four different two 



particle energies. We notice that the set of coupled Eqs. 
((39)) for the evolution of the populations and of the spin 
accumulation has a similar structure to that reported 
jj ^29^3i,39 fpj. single level quantum dot, a metallic island 
and a single-walled carbon nanotube, respectively. Some 
prefactors and the argument of the principal part terms, 
however, are DD specific. In particular, as i n^9,39 ^ 
clearly identify a spin precession term originating from 
the combined action of the reflection at the interface and 
the interaction. The associated effective exchange split- 
ting is 7i3i, with 7 = —.9 /is being the gyromagnetic ra- 
tio, and 



(42) 



being the corresponding effective exchange field. We fo- 
cus now on the stationary limit. In the absence of the 
precession term the spin accumulation has only a S^J''^ 
component since, due to our particular choice of the spin 
quantization axis, — holds. The exchange field tilts 
the accumulated spin out of the magnetizations' plane 
and gives rise to a nonzero component proportional 

to Bi and 4^^- 

To get further insight in the spin-dynamics we observe 
that, since we are looking at the low voltage regime, we 
can linearize the Fermi function fa in the bias voltage, 
i.e., 

= (l + e''(«+^^=))-i « /(0(l-/(-Oe/3X4). (43) 
Introducing the polarization of the contacts 



DauiO+DalAO 



we can express the F 



factors as 



(0 



where — D^^^ + Dai^- It is also sufficient for our 
calculations to regard the density of states as a constant 
quantity, Da{^) = D^- Consequently the polarization is 
also constant, PaiO — Pa- Finally, we focus in the fol- 
lowing on the symmetric case where both leads have the 



same properties, which in particular means that tunnel- 
ing elements, polarizations, density of states and reflec- 
tion amplitude are equal: 



t, Pl = P2 



Di — D2 :— D, i?i(jQ_o-Q — i?2cr0-CTo := R- (44) 

Upon introducing the lincwidth F = the conduc- 

tance G12 — 112/Vbias for the resonant regime iV = 1 ^ 
N ^2 reads 



r (P,\ - ^-2>?f^2 /(M2)/(-M2) 

p" sin^d) 



(45) 



\^ l + [i3l//(/i2)2rfc2]2cos2(f)y 

Similarly we find for an arbitrary resonance (i = 0, 1, 2, 3) 

G,.+r(e) = \e^PW + (46) 
2 l + /((-l)>,+i) 

P^sin^(l) \ 

1^ l + [B,+i//((-l)^+V.+i)2r|(z + lMt|z)|2]2cos2(f)j 

where |(« + l|(i^^|i)| is a shortcut notation for the non 
vanishing matrix elements -I- l|d]jQ|i?|*'''z)| calcu- 

lated in the tables of Appendix II. It holds KIMIqIO)] = 

|(4|4q|3)| = I/V2, and |(2|4q|1)| = |(3|40|2)| = k+. 
Moreover, we gathered together the principal part contri- 
butions and the ones coming from the refiection Hamil- 
tonian in the effective magnetic fields 



Bo 



B, = B4:^-^|iJX(A*4,i?f -{£^2})^ 
The latter are defined in terms of the function 



{E2}) 



[P+^^{^i4) + Rc.u] 



~^^«n(^3o - £2") - k^^P^ni^io - £2'")- 

Moreover, a closer look to Eq. shows that its angu- 
lar dependence is strongly coupled to the square of the ra- 
tio {'-fBi)/{hTf{{—iy^i)), which is the effective exchange 
splitting, rescaled by the coupling and the Fermi func- 
tion. The ratio occurs in the denominators, and its value 
depends on the gate voltage. As the change of Bi un- 
der variation of the gate voltage is comparatively small, 
the factor dominating the gate voltage evolution is the 
Fermi function. This accounts for the population of the 
dot: only if a nonzero spin is present (i.e. odd filling: one 
or three electrons), the effective magnetic field can have 
an infiuence. That is why correspondingly the renormal- 
ized effective exchange splitting vanishes for even fillings. 
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FIG. 3: (Color online) Gate voltage dependence of the factors 
Bf / f^r^ entering the conductance formula (|46p . Notice the 
mirror symmetry of the 0^1 with the 3 4 curve and of 
the 1^2 with the 2^3 one. 



namely below the <-^- 1 and 2^3, respectively above 
the 1^2 and 3 4 resonances. The curves belonging to 
the resonances involving the half filling do not go imme- 
diately to zero but show a more complex behavior with 
some small intermediate peaks due to the influence of 
the various excited states present for a two-electron pop- 
ulation of the dot. This can nicely be seen from figure 
[3] (remember that Vgate —^), where the four different 
factors (7i?i)^/(/((— l)*/ii)r)^ are plotted. As we expect 
Goi(0 and G'34(^), respectively Gi2(0 and G23{£,) are 
mirror symmetric with respect to each other when the 
gate voltage is varied. This in turn reflects the electron- 
hole symmetry of the DD Hamiltonian. The parameters 
of the figures are chosen to be (6 < 0) 



U 



4-10"2|6|, hT ^ 
6\bl V = 1.6\bl 



4- 10-3|5|, 



(47) 



and p = 0.8, B = 0.05Z?. As expected, the peaks are mir- 
ror symmetric with respect to the half-filling gate voltage. 
Notice also the occurrence of different peak heights, both 
in the parallel as well as in the antiparallel case. For both 
polarizations the principal-part-terms entering Eq. (I46p 
vanish, the spin accumulation is entirely in the magneti- 
zation plane and the peak ratio is solely determined by 
the ratio of the groundstate overlaps 2/fc^. For polariza- 
tion angles ^ 0, tt the ratio is also determined by the 
non-trivial angular and voltage dependence of the effec- 
tive exchange fields. Finally, as expected from the con- 
ductance formulas (^5]) . the conductance is suppressed 
in the antiparallel compared to the parallel case. The 
four conductance peaks are plotted as a function of the 
gate voltage in Fig. [4] for the polarization angles 8 = 
and 8 = TT, top and bottom figures, respectively. These 
features of the conductance are nicely captured by the 
color plot of Fig. O where numerical results for the con- 
ductance plotted as a function of gate voltage and polar- 
ization angle are shown. The conductance suppression 
nearby 8 = tt is clearly seen. 

In the following we analyze in detail the single res- 
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FIG. 4: (Color online) Conductance at low bias for the par- 
allel case = (top) and antiparallel case Q = tt (bottom). 
Notice the different peak heights and the mirror symmetry 
with respect to the half filling value ^ = 0. The conductance 
in the antiparallel configuration is always smaller than in the 
parallel one. 



G[eyh] 




FIG. 5; (Color online) Conductance as a function of the po- 
larization angle and of the gate voltage. The minimal con- 
ductance peaks occur as expected at = tt 



onance transitions. Due to the mirror symmetry it is 
convenient to investigate together the resonances N = 
0^N=1, N = 3^N = AsindN=l^N = 2, 
N — 2 ^ N — 3. We use the convention that, for a fixed 
resonance, the parameter ^ = when /iyv+i = 0. 



1. Resonant regimes N = ^ N = 1 and N = 3 ^ N = 4 

The expected mirror symmetry of Gqi and G34 is 
shown in Figure[Bl where the conductance peaks are plot- 
ted for different polarization angles 8 of the contacts. 
Notice that the analytical expressions (continuous 
lines) perfectly match the results obtained from a nu- 
merical integration of the master equation (|2ip with the 
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FIG. 6: (Color online) The conductance Goi(C) (upper fig- 
ure) resp. G34(^)) (lower figure) vs gate voltage for different 
polarization angles. The mirror symmetry of the conductance 
peaks for the 1 and 3^4 transitions is clearly observed. 
Notice the excellent agreement between the prediction of the 
analytical formula Eq. (|46|) (continuous lines) and the results 
of a numerical integration of Eqs. (|21|l with (|25p (symbols). 
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FIG. 7: (Color online) Tunneling magnetoresistance (TMR) 
for the transitions 0^1 (upper figure) and 3^4 (lower 
figure) vs gate voltage. The TMR is always positive and is 
independent of gate voltage for coUinear lead magnetizations 
only, = and Q = -k. 



current formula ([25)1 . We also can see that the maxima 
of the conductance decrease with growing up to tt. It 
can be shown that the peaks for = and = tt he at 
the same value of ^, because the effective fields Bi exactly 
vanish due to trigonometrical prefactors. In other words, 
virtual processes captured in the effective fields Bi do not 
play a role in the colhnear case. For noncoUinear config- 
urations, however, the peak maxima are shifted towards 
the gate voltages where an odd population of the dot 
dominates, because there the effective exchange field can 
act on the accumulating spin and makes it precess, which 
eases tunneling out. These findings are in agreement with 
results obtained for a single- level quantum dot'^^ , a metal- 
lic island^- and carbon nanotubes"^^. To quantify the rel- 
ative magnitude of the current for a given polarization 
angle O with respect to the case 9 = we introduce the 
angle-dependent tunneling magnetoresistance (TMR) as 



TMi?^w+i(e,0 = 1 



G 



JV,JV+1 



G 



(0,6 



For the transition ^ 1 it reads 
TMRm 



l + [i?2//2(-^i)r2]cos2(f) 



(48) 



Hence, the TMR vanishes for 8 = and takes the con- 
stant value 
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FIG. 8; (Color online) Goi(e)/Goi(0) (upper figure) resp. 
G34(0)/G34(O) (lower figure) vs polarization angle. For all 
the three chosen values of the gate voltage the curve display an 
absolute minimum at O = tt. Notice the overall agreement of 
the analytical predictions (|46[) given by the continuous curves 
with outcomes of a numerical solution of the master equation, 
Eq. ([21} , together with (symbols) . 
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FIG. 9: (Color online) Gi2(e)/Gi2(0) (upper figure) resp. 
G23(0)/G23(O) (lower figure) vs polarization angle. Notice 
the overall agreement of tfie analytical predictions (|45|) and 
(|46|) . continuous lines, with the data (symbols) coming from 
a numerical solutions of the equations for the reduced density 
matrix. 



at = TT. For the remaining polarization angles, & ^ 
and 6 7^ TT, the TMR is gate-voltage dependent and pos- 
itive. The behavior of the TMR as a function of the 
gate voltage is shown in figure [T] To understand the gate 
voltage dependence of the TMR at noncoUinear angles we 
have to remember that the dot is depleted with raising 
^. For the transition <-> 1, this means that at positive 
^, the dot is predominantly empty, so that an electron 
which enters the dot also fast leaves it. In this situation 
the TMR is finite and its value depends in a complicated 
way on the amplitude of the exchange field. At negative 
5, the DD is predominantly occupied with an electron 
which can now interact with the exchange field, which 
makes the spin precess and thus eases tunneling out of 
the dot. Consequently, G'ArAr+i(9, <^) « Gnn+i{0,S,) and 
the TMR vanishes. Finally, figure [8] illustrates the angu- 
lar dependence of the normalized conductance for three 
different values of the gate voltage. We detect a com- 
mon absolute minimum for the conductance at = tt, 
i.e., transport is weakened in the antiparallel case. The 
width of the curves is dependent on the rcnormalizcd ef- 
fective exchange {'^Bi) / {T f). The larger its value the 
narrower get the curves, because the spin precession can 
equilibrate the accumulated spin for all angles but = tt. 
Notice again the equivalence of the curves belonging to 
^ = ±2|&| for the 1 <-> 2 resonance to the curves with 
^ = t2|&| for the 3 4 resonance. 
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FIG. 10; (Color online) Differential conductance ^ for the 
parallel B = (top), perpendicular O = 7r/2 (middle) and 
antiparallel B = tt (bottom) configurations. The two half di- 
amonds and three diamonds regions correspond to bias and 
gate voltage values where transport is Coulomb blocked. The 
excitation lines, where excited states contribute to resonant 
transport, are clearly visible in all of the three cases. How- 
ever, a negative differential conductance is observed in the 
perpendicular case, while some excitation lines are absent in 
the antiparallel configuration. 



2. Resonant regimes N — 1 ^ N — 2 and N 



N : 



For the resonant transitions 1^2 and 2 3 qual- 
itatively analogous results as for the 0^1 and 3^4 
transitions are found. Thus, exemplarily we only show 
the angular dependence of the normalized conductance 
in Fig. [21 showing the expected absolute conductance 
minimum at = 0. 



V. NONLINEAR TRANSPORT 

In this section we present the numerical results, de- 
duced from the general master equation (j2ip combined 
with the current formula (1251). We show the differen- 



tial conductance ^(^, Vbias) for the three distinct angles 
Q — 0, Q — ^ and Q — n, see figure [TOl top, middle and 
bottom, respectively. The results confirm the electron- 
hole-symmetry and the symmetry upon bias voltage in- 
version I{^,VMas) = -Vbias)- In aU of the three 
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FIG. 11: (Color online) Current (top) and triplet occupation 
(bottom) for the two coUinear {O = 0,0 = tt) cases and the 
perpendicular case (O = |-) at a fixed gate voltage ^ = 4|6|. 
Notice the occurrence of a pronounced negative differential 
conductance feature for perpendicular polarization O — -k 12. 



cases we can nicely see the expected three closed and 
the two half-open diamonds, where the current is blocked 
and the electronic number of the double-dot-system stays 
constant. At higher bias voltages the contribution of ex- 
cited states is manifested in the appearance of several 
excitation lines. One clearly sees that transition lines 
present in the parallel case are absent in the antiparallel 
case. Moreover, in the case of noncoUinear polarization, 
Q = tt/2, negative differential conductance (NDC) is ob- 
served. 

In the following, we want not only to explain the origin 
of these two features, but alongside also give another ex- 
ample for spin-blockade effects, which play a decisive role 
in the DD physics. 

As a starting point, we plot in figure [TT] (top) the cur- 
rent through the system for the three different angles 
Q = |0, ^, tt} at a fixed gate voltage ^ = 4|6| and posi- 
tive bias voltages. We recognize, that for eVsias < 2.4|6| 
the current is Coulomb-blocked in all of the three cases. 
In this configuration exactly one electron stays in the 
double-dot. From about eVsias > 2.4|6| the channel 
where the groundstate energies fii and /i2 are degenerate 
opens (I leer) <-> |2) transition) and current begins to flow. 
With increasing bias more and more transport chan- 
nels become energetically favorable. In particular, for 
all the polarization angles Q we observe two consecutive 
steps corresponding to the transitions |0) ^ I lee) and 
|lecr) |2'). The latter, occurring at about eVbias — 4|6|, 
involves the excited two-particle triplet states \2'{Sz))- 
The next excitation step, indicated with a circle in Fig. 
[TT] (top), belongs to the transition \loa) ^ |2"). The 
associated line is missing for the antiparallel configura- 
tion, as well as the lines corresponding to \loa) ^ |2"'); 
|lecr) ^ |2"); |2) ^ \3oa); \lea) ^ |2"'). Crucially, 



in all of these transitions a two-particle state with total 
spin zero is involved. In order to explain the absence 
of these lines, let us e.g. focus on the first missing step 
corresponding to the \loa) <-*■ |2") resonance. In the par- 
allel case (say both contacts polarized spin-up) there is 
always an open channel corresponding to the situation in 
which the spin in the DD is antiparallel to that in the 
leads (i.e. |lo— )). In the antiparallel case (say source 
polarized spin-up, drain polarized spin-down) originally 
a spin-down might be present in the dot. An electron 
which enters the DD from the source must then be spin- 
up (in order to form the state |2")), but as the drain is 
down-polarized, it will be the spin-down electron which 
leaves the DD, which corresponds to a spin flip. Now 
the presence of a spin-up electron in the DD prevents a 
majority (another spin-up) electron from the source to 
enter the DD, such that we end up in a blocking state. 
The transition is hence forbidden. 

A similar, yet different spin-blockade effect determines 
the occupation probabilities for the triplet state. Fig. 
[TT] (bottom). Naturally, for all angles the probability 
to be in the triplet state increases above the resonance 
at eVbias — 4|6|, but interestingly, such probability is 
largest in the antiparallel case. This is due to the fact 
that a majority spin in the parallel configuration (spin- 
up) can be easily transmitted through the DD via the 
triplet states |2'(1)) or |2'(0)). In the antiparallel case, 
however, a blocking state establishes (say again source 
polarized spin- up, drain polarized spin-down). Let ini- 
tially a spin-down electron be present on the DD. From 
the source electrode, most likely a majority electron po- 
larized spin-up electron will enter the dot. Now, just as 
in the previous case, the consecutive tunneling event will 
cause a spin flip in the DD, because the spin-down elec- 
tron (majority electron of the drain) will leave the dot. 
So the DD is finally in a spin-up state, and once the next 
majority spin-up electron from the source enters, the DD 
ends up in the triplet state |2'(-|-1)) and will remain there 
for a long time due to the fact that the majority spins 
in the drain are down-polarized. Hence the triplet state 
|2'(-|-1)) acts as a trapping state. 

Notice that the two distinct spin-blockade effects are 
different from the Pauli spin-blockade discussed in the 
DD literature^^'"'^'^^"''". Moreover, the second effect, re- 
lying on the existence of degenerate triplet states, is also 
different from the spin-blockade found in Refi^ for a sin- 
gle level quantum dot. 

Finally, let us turn to the negative differential con- 
ductance, which occurs for noncollinearly polarized leads 
(see the dashed blue lines in figure [TT]) , and which we 
find to become more evident for higher polarizations (not 
shown). Neglecting the exchange field, we would just ex- 
pect the magnitude of the current for the noncoUinear 
polarizations to lie somewhere in between the values for 
the parallel and the antiparallel current, because the non- 
coUinear polarization could in principle be rewritten as 
a linear combination of the parallel and the antiparal- 
lel configuration. Now the effect of the exchange is to 
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cause precession and therewith equilibration of the accu- 
mulating spin, which corresponds to shifting the balance 
in favor of the parallel configuration, i.e. enhancing the 
current. The decisive point is that the exchange field 
is not only gate, but also bias voltage dependent and 
reaches a minimum around eVsias ~ 8\b\. This explains 
the decreasing of the current up to this point. Afterwards 
the influence of the spin precession regains weight. The 
same consideration applies for the other NDC regions ob- 
served in Fig. [TUl e.g. in the gate voltage region « 2\b\ 
involving the = iV = 1 transition, as described in 
Ref^. 



Wn 



where 



-1 E [F^i(Pn)W^, 

a—s,d 



MiT/i = Ml T Eze 



(50) 



(51) 



and F^^^^{E) = DaaQf^{E). For the collinear case is 
DaaQ = Da±^ if <Jq =] / i in the parallel case. On the 



other hand D 



S(JQ 



/ I, in the antiparallel case. 



Ds^^ and D^aQ 



D 



if (Tr. 



VI. THE EFFECTS OF AN EXTERNAL 
MAGNETIC FIELD 

In this section we wish to discuss the qualitative 
changes brought by an external magnetic field applied to 
the DD. Specifically, the magnetic field is assumed to be 
parallel to the magnetization direction of the drain. For 
simplicity we focus on the experimental standard case of 
parallel and antiparallel lead polarization and of low bias 
voltages. Then, the magnetic field causes an energy shift 
^Ezeeman depending on whether the electron spin is par- 
allel or antiparallel, respectively, to it. For collinear po- 
larization angles the principal part contributions vanish, 
and the equations for the RDM are easily obtained. We 
report exemplarily results for the transitions <-> 1 and 
1 2. Let us then consider the parameter regime nearby 




FIG. 12: (Color online) Conductance vs gate voltage for par- 
allel (continuous line) and antiparallel (dashed lines) contact 
configuration and Zeeman splitting Ezeeman ~ 0.05|6|. The 
magnetic field breaks the mirror symmetry with respect to 
the gate voltage in the parallel configuration. 



the 0^1 resonance, and setup a system of three equa- 
tions with three unknown variables W^OiW^iT ^^'^ ^il- 
The first equation corresponds to the normalization con- 
dition Wii + W11+W0 — 1. The remaining equations are 
the equations of motion for W^i|/x which can be written 
as 



lie I 



^aT(MiT)W^iT 



IT 



(49) 
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FIG. 13: (Color online) Tunneling magnetoresistance (TMR) 
vs gate voltage in the presence of an external magnetic field. 
In contrast to the zero field case, the TMR can become neg- 
ative in the vicinity of the 2 3 and 3 4 resonances. 



Upon considering symmetric contacts (^i — t2 — t, 
Di = D2 = D) we find in the parallel case 

Goi(e = 0) = ^/(-A*iT)/(-mi) (52) 

^ P (/(Mit) - /(Mu)) + /(Mit) + /(Mu) 

/(-Mit)/(Mu) + /(-Mit)/(-Mu) + /(MiT)/(-Aiu) ' 

For the antiparallel case we obtain 



Goi(e = 7r)=Goi(e = 0) 

i-p^(/(m) + /(Mu)) 

p[/(Mit)-/(mu)]+/(Mit) + /(mu)' 
Analogously we find for the 1 <-> 2 transition 
Fe^fci 



(53) 



Gi2(e = o) 



2/3 



-/(M2T)/(M2i) 



(54) 



p(/(-M2x) - /(-M2t)) + f{-m) + /(-M2i) 

/(-/i2T)/(M2i) + /(-M2T)/(-M2i) + /(M2T)/(-M2i) ' 



Gi2(e-^)-Gi2(e = o) (55) 

^ l-P^(/(-M2T)+/(-M2i)) 

p > [/(-/"2i) - /(-M2t)] + /(-M2t) + f{-m) ' 

The remaining resonances are calculated analogously. 
Figure [T^ shows the four conductance resonances for the 
parallel and antiparallel configurations. Strikingly, the 
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applied magnetic field breaks the symmetry between the 
tunneling regimes <-^- 1 and 3 <-> 4, as well as between 
the resonances 1^2 and 2 ^ 3 in case of parallel contact 
polarizations. The reason for this behavior is the follow- 
ing: In the low bias regime transitions between ground- 
states dominate transport. In particular, the magnetic 
field removes the spin degeneracy of the states {lea) and 
of the states |3o(t), such that states with spin aligned 
to the external magnetic field are energetically favored. 
Therefore, the transport electron in the tunneling regime 
<-^- 1 is a majority spin carrier. For the case 3 <-> 4, how- 
ever, two of the three electrons of the groundstatc |3o+) 
have spin-up, such that the fourth electron which can be 
added to the DD has to be a minority spin carrier. There- 
fore, the conductance gets diminished with respect to the 
1 transition, and the mirror symmetry present in the 
zero field case is broken. Analogously the broken sym- 
metry in the case of the transitions 1 <~+ 2 and 2 <-> 3 can 
be understood. Correspondingly, the TMR can become 
negative for values of the gate voltages around the 2 3 
and 3 4 resonances. We observe that a negative TMR 
has recently been predicted in^ for the case of a single 
impurity Anderson model with orbital and spin degener- 
acy. In that work a negative TMR arises due to the as- 
sumption that multiple reflections at the interface cause 
spin-dependent energy shifts. In our approach, however, 
where the contribution from the reflection Hamiltonian is 
treated to the lowest order, see (fT9)) . such spin-dependent 
energies are originating from the magnetic field-induced 
Zeeman splitting. 

VII. CONCLUSIONS 

In summary we have evaluated linear and nonlinear 
transport through a double quantum dot (DD) cou- 
pled to polarized leads with arbitrary polarization di- 
rections. Due to strong Coulomb interactions the DD 
operates as a single-electron transistor, a F-SET, at low 
enough temperatures. A detailed analysis of the current- 
voltage characteristics of the DD, and comparison with 
results of previous studies on other F-SET systems with 
noncoUinear (a single level quantum dot^, a metallic 
island^^, a carbon nanotubc'^^), bring us to the identi- 
fication of universal behaviors of a F-SET, i.e., a be- 
havior shared by any of those F-SETs, independent of 
the specific kind of conductor is considered as the central 
system, as well as system-specific features. 
Universal is the presence of an interfacial exchange field 
together with an interaction-induced one, present for 



noncoUinear polarization only. These exchange fields 
cause a precession of the accumulated spin on the dot and 
therewith ease the tunneling out. This effect has various 
implications. It determines e.g. the gate and angular 
dependence, as well as the height of single conductance 
peaks and can yield negative differential conductance fea- 
tures. Another universal feature is the occurrence of a 
negative tunneling magnetoresistance - even in the weak 
tunneling limit - if an external magnetic field is applied. 
Specific to the DD system are the following features: In 
the low bias regime, the problem can be solved analyt- 
ically and for the tunneling regimes <-> 1 and 1 <-> 2, 
respectively 2 <-> 3 and 3 <-*■ 4, the system behaves equiv- 
alent to a single level quantum dot, where the Coulomb 
blockade peaks are found to be mirror-symmetric with 
respect to the charge neutrality point. This mirror sym- 
metry reflects the electron-hole-symmetries of a system 
and is therefore typical for the DD, as well as the ra- 
tio of the peak heights. An external magnetic field lifts 
this symmetry and can cause a negative tunneling magne- 
toresistance. In the nonlinear bias regime, the presence of 
various excited states gives rise to interesting DD specific 
features. For example, a suppression of several excitation 
lines for an antiparallel lead configuration originates from 
a spin-blockade effect. It occurs because a trapping state 
is formed whenever a transition involves a two-electron 
state with total spin zero. A second spin-blockade ef- 
fect we described involves the two-electron triplet state. 
The common mechanism of these two spin-blockades is 
the following: in both cases, a tunneling event can only 
occur if initially the dot is populated with an unpaired 
electron possessing the majority spin of the drain. The 
second step is that a majority electron of the source will 
enter, forming a spin-zero state. Then the first electron 
can leave the dot, causing a spin-flip. If the triplet state 
is involved, the dot will be left in a trapping state once 
a second majority electron from the source enters. Oth- 
erwise, we are directly in a blocking state. Finally, for 
noncoUinear lead polarizations, negative differential con- 
ductance can be observed. 

All in all, due to their universality and the multiplicity of 
their properties, F-SET based on DD systems seem good 
candidates for future magneto-electronic devices. 
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APPENDIX A: MATRIX ELEMENTS OF THE 
DOT OPERATORS 

The following tables show all the possible matrix- 
elements {N - l|da<70 \N) and (A^MLq - 1) with a = 
{1,2} and ctq — {T)i} which occur in the master equa- 
tions (|Bip and (|B2p . Notice that we only need to il- 
lustrate cither the matrix elements {N — l\daaQ\^) or 
{N\dl^^^\N — 1), because they are complex conjugated 
to each other. 
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TABLE II: Matrix elements for the transition iV = ^ iV = 1 
induced by the operators dot and d^i, a = 1,2. 
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TABLE III: Matrix elements for the transition iV = 1 ^ iV = 
2 induced by and dj^^, a = 1,2. The notation |2'(sz)), 
with = 0, ±1, specifies which one of the triplet elements is 
addressed. 
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TABLE IV: Matrix elements for the transition iV = 2 ^ iV = 
3 governed by dat and d^i, a = 1,2. 
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TABLE V: Matrix elements for the transition A = 3 ^ A = 4 
induced by the operators d^. and dj^, , a = 1, 2. 
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APPENDIX B: THE MASTER EQUATION FOR 
THE RDM IN THE LINEAR REGIME 

We report here explicitly the coupled equations of mo- 
tion for the elements and pmrt^\t) of the RDM 
to be solved in the low bias regime. They are obtained 
from the generalized master equation (|2ip upon observ- 
ing that i) in the linear regime terms that couple states 



with particle numbers unlike N and N + I can be ne- 
glected; ii) we can reduce the sum over h and h' and 
over I, /' only to energy ground states. In the remain- 
ing not energy conserving terms the sum has to go also 

over excited states. With ^n+i ■— e'^^+i — e'^^ being 
the chemical potential we finally arrive at the two master 
equations 



(Bl) 



4 E E E - E 
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Notice that we kept the sums over excited states I, h in transitions. 
(|B1I) and l,hm (|B2p which are responsible for the virtual 
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